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In this chapter we discuss methods to calculate the entanglement of a system using density-functional theory. 
We firstly introduce density-functional theory and the local-density approximation (LDA). We then discuss 
the concept of the 'interacting LDA system'. This is characterised by an interacting many-body Hamiltonian 
which reproduces, uniquely and exactly, the ground state density obtained from the single-particle Kohn-Sham 
equations of density-functional theory when the local-density approximation is used. We motivate why this idea 
can be useful for appraising the local-density approximation in many-body physics particularly with regards to 
entanglement and related quantum information applications. Using an iterative scheme, we find the Hamiltonian 
' characterising the interacting LDA system in relation to the test systems of Hooke's atom and helium-like atoms. 

, The interacting LDA system ground state wavefunction is then used to calculate the spatial entanglement and 

. the results are compared and contrasted with the exact entanglement for the two test systems. For Hooke's atom 

' we also compare the entanglement to our previous estimates of an LDA entanglement. These were obtained 

using a combination of evolutionary algorithm and gradient descent, and using an LDA-based perturbative 
approach. We finally discuss if the position-space information entropy of the density — which can be obtained 
directly from the system density and hence easily from density-functional theory methods — can be considered 
as a proxy measure for the spatial entanglement for the test systems. 

I. INTRODUCTION 

■ The mapping of an interacting many-body system to a collection of effective single particle equations is a ubiquitous and 
powerful approach in theoretical physics. This is perhaps most notably embodied in density-functional theory (DFT) where 
r a non-interacting many -body system, the Kohn-Sham system is proven to have exactly the same ground state density as 
■ O the original interacting system. By solving a set of single-particle Kohn-Sham (KS) equations, the ground state density can be 
calculated and from this in principle all the ground-state properties of the interacting many -body system |2]. The KS equations 
would give the exact density due to the inclusion in their effective single-particle potential of the so-called "exchange-correlation" 
potential (Wxc)- In practice is an unknown functional of the density so approximations for Wxc have been designed. The most 
popular and often effective method is the local-density approximation (LDA). DFT has been so far one of the most powerful 
methods to calculate many-body properties of complex systems, and its efficiency lies in the observation that modelling many- 
particle systems exactly is generally computationally intractable as, for example, one is dealing with a many-body wavefunction 
whose storage requires excessive amounts of memory. Solving a system of single-particle equations or dealing with the density, 
which is a function of only one variable, represents indeed a computationally attractive alternative. However we do not know in 
P\| general how to express some ground state properties, such as the entanglement, in terms of the ground state particle density, and 
this implies that common approximation schemes used within DFT do not have a straightforward extension to the calculations of 
these properties. Here we present a scheme to calculate the entanglement associated with the LDA approximation and investigate 
how well it compares with the exact entanglement. Our test-bed systems are helium-like atoms and Hooke's atoms. This 
T— I . information could be useful when considering which approximation to invoke to model the suitability of systems as a component 
' of a quantum information processing device. As the entanglement is considered as one of the main resources for these devices, 
the knowledge of how much entanglement could be created in a certain system appears essential for quantum information 
purposes. In addition comparison of the entanglement associated with a DFT approximation to the exact entanglement offers 
another way to appraise the accuracy of the approximation. 

In this chapter we firstly introduce DFT and the LDA. We explain the concept of the 'interacting LDA system' (i-LDA) and 
how this may be used to calculate the entanglement corresponding to the use of the LDA. We then introduce our test-bed systems, 
Hooke's atom and the helium-like atoms. Next we consider how these systems can be solved using numerically 'exact' methods 
and also approximately using the LDA within DFT. The exact and i-LDA spatial entanglement results are then displayed and 
compared. We contrast these results with the position space information entropy which can easily be expressed as a functional 
of the density — and hence calculated using traditional DFT schemes — and investigate its usefulness as a proxy measure for the 
entanglement. Finally, we compare the i-LDA entanglement to our previous best approximations to an LDA entanglement, found 
using a combination of evolutionary algorithm and gradient descent, and an LDA-based perturbative approach | 
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II. DENSITY-FUNCTIONAL THEORY 



Directly solving the many-body Schrodinger equation is a formidable, arduous task and doing so accurately for systems of 
more than a few electrons remains outside the realms of contemporary computing. This can be partly illuminated by considering 
the wavefunction on a mesh of 1000"^ points for one electron; for the same accuracy with N electrons one would require 1000'^^ 
points, clearly an excessive amount of memory as N becomes large. Yet nature, in a sense performs similar calculations all 
the time, but these are too complex to be reproduced by 'classical' computation. Quantum computing offers some hope in the 
simulation of quantum systems (see for example Ref. |i4J and Ref. [5J). However, unless quantum computers of a substantial 
number of qubits are realised, simulating large systems may not be feasible. There is though a powerful and important method 
that allows classical computing to efficiently model many-body electron systems through the density: density-functional theory. 
The density depends upon only one co-ordinate (so only 1000"^ mesh points in our example) therefore by using this, rather 
than the many-body wavefunction, solving many-body systems can be made computationally more tractable. In its simplest 
formulation density-functional theory expresses observables of the ground-state of an interacting electron system as functionals 
of the ground-state density. 



with a non-degenerate ground state, the ground state density uniquely determines the ground state many-body wavefunction 
and therefore the external potential. This is known as the Hohenberg-Kohn theorem |2]. The usual proof shows how the 
density uniquely determines the wavefunction |6]. Assume there exist two Hamiltonians Hi and H2 with different ground-state 
wavefunctions, and [^'2) respectively, that give rise to the same density. Then 




(1) 



and is exact in principle due to a theorem by Hohenberg and Kohn which we present next. 



A. The Hohenberg-Kohn Theorem 



In an electron system with Hamiltonian 




(2) 



(^-lli/zl^l) = (*l|if2 
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Then by the variational principle 




(4) 



with E2 the ground state related to H2. 

This can be repeated but with 1 and 2 exchanged to give 



{^2\Hi 1*2) = (*2|-ffl 
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(wi,ext(r) - V2,ext{r)) n{r)dr + E2. (5) 



By the variational principle 



J {vi,ext{r) - V2^ext{r)) n{r)dr + E2 > Ex. 



(6) 



Adding the inequalities gives 



Ei+E2> E1+E2. (7) 



This inequality is a fallacy so the initial assumption must have been incorrect: there do not exist two ground state wavefunctions 
that give rise to the same density. Hence by reductio ad absurdum the density uniquely determines the wavefunction. 
As we may express the external potential as 

i 

then the density by determining the wavefunction also determines the external potential (up to an additive constant), and therefore 
all of the properties of the system. 

Although the Hohenberg-Kohn theorem proves the existence of the one-to-one relationship between ground state densities 
and wavefunctions it does not suggest how to construct the density nor find the ground-state properties of the system from the 
density. To achieve this the Kohn-Sham equations iJl are used. 

B. The Kohn-Sham equations 

The Kohn-Sham (KS) equations allow the ground-state density to be found efficiently and the system energy to be expressed 
using the density. The KS equations are derived by considering the system energy as a functional of the density ^^[n], this is 
minimised subject to a constraint incorporating N, the number of electrons, 

j n{r)dr = N. (9) 
Using Lagrange multipliers with the functional derivative 

5F[n{y)] ^ ^.^^ F[n{y)+e5{y-x)]-F[n{y)] ^^^^ 
5n{x) e^o e 



5E\ri\ 

^-M = 0. (11) 

on[r) 

As the expression of the kinetic energy as a functional of the density for an interacting electron system is unknown, Kohn and 
Sham re-wrote the exact energy functional in terms of the kinetic energy of non-interacting electrons TNi[n], the Hartree energy 
U\n\, and the potential energy from the electron density. This requires the introduction of a functional correctional term, the 
exchange-correlation energy Exc[n]- The expression for E[n\ then becomes 

E[n] = TNi[n] + U[n] + / v,.,t{r)n{r)dr + E,,[n] (12) 
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where 



Using Lagrange multipliers now gives 

5TNi[n] 



+ Veffir) - (1 = 0. (14) 



Sn{r) 
Here 

'"eff{r) = Viiir; [n]) + Vext{r) + v^^{r; [n]) (15) 
where they define the exchange-correlation potential 



and the Hartree potential 



Vy.c{r\[n\) = —-— (16) 
ori(r) 



^H(r;W)= / ^^dr'. (17) 
' r — r'\ 



Now Eq.[T4]can be read as a non-interacting system — the KS system — with a potential v^ff- So to find the density and hence 
the energy the KS equations must be solved self-consistently: 

-iv2 + «,/^(r)^ 0,(r) = e,0,(r) (18) 
and 

N 



n(r) = ^|0,(r)|2. (19) 

i=l 

As 

TnM =Y1J '^n^) (-^^''^*(^)) ^20) 

then the energy E[n] may be calculated from the density using E^c- 

This scheme would be exact but the general functional form of E^c is unknown hence it has to be approximated. One of 
the simplest yet often effective approximations is using the fact that, for systems with slowly varying density, the exchange- 
correlation density can be locally approximated by that of a uniform electron gas. This is known as the local density approxima- 
tion. 

C. The local-density approximation 

The local-density approximation for the exchange-correlation energy depends only upon the density as opposed to generalised 
gradient approximations which include first and higher derivatives. It may be written using the form 



E^^ = j dre^cHr)). (21) 



The exchange-correlation energy density may be separated into its exchange and correlation parts, Cxc = + Cc- Then, for a 
three-dimensional homogeneous electron gas, Cx is known exactly \]^ 

3 /3^^/^ 



ex = --[-] n^, (22) 
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while Be has to be approximated. We use Cc as calculated numerically by Perdew-Wang fs*] by fitting the parameters A, a, 
Pi.. to the results of Monte-Carlo simulations of the homogeneous electron gas |9]. For zero spin polarisation this has the 
form 

ec = ~2An{l + ars)\og\l + ^ ^- ^ ), (23) 

where is the radius of the sphere whose volume is equal to the volume per electron in a homogeneous electron gas, 

1 
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47rn(r) 



(24) 



The above approximation is exact for the homogeneous electron gas so would be expected to be more accurate for systems close 
to this limit. 



III. THE INTERACTING LDA SYSTEM 



We wish to compare the entanglement calculated using a local-density approximation to the exact entanglement, and although 
in theory the density determines all the ground-state properties of the system it is unknown how to use the density to represent 
certain observables, including the entanglement. One method to calculate the entanglement within LDA is to construct that 
many-body interacting system whose ground state wavefunction reproduces exactly the LDA density. We named this 'the 
interacting LDA system' (i-LDA) IH. By the Hohenberg-Kohn theorem |01 this system is unique. 

Given a test-bed interacting many-body system for which ground-state properties, such as the entanglement, can be calculated 
exactly, by comparing their values with the ones obtained using the corresponding i-LDA wavefunction, we can infer how precise 
the LDA can be as an approximation to those quantities. Notice that the external potentials characterising the initial many-body 
system and the i-LDA system are different, and that the LDA density is at the same time an approximation to the density of the 
test-bed system and the exact ground state density of the i-LDA system. 

Such a route to assess system properties within the LDA is computationally intensive, but, for small test systems, allows a 
way to investigate the accuracy of the LDA for quantities for which functional expressions in terms of the density are currently 
unknown. In addition this procedure offers a novel way to appraise the LDA by considering the differences between the exact and 
i-LDA external potentials, thereby furnishing researchers with spatially resolved information about an approximation's accuracy. 
A similar procedure to the one described can be applied to other approximations within DFT. 

We introduced a method to approximate the i-LDA wavefunction in Ref. iH, then in Ref. ifToll we developed an iterative 
scheme to find the Hamiltonian of the i-LDA system 'exactly' (within numerical accuracy). Then the 'exact' i-LDA wavefunction 
could then be derived. This is the scheme we use in this chapter to allow the entanglement corresponding to the LDA to be 
ascertained. 



rV. THE TEST-BED SYSTEMS 



The helium atom comprises two electrons held by a nuclear charge arising from two protons in the nucleus. By using the Born- 
Oppenheimer approximation the motion of the nucleus may be ignored. The system generalises to helium-like atoms by consid- 
ering a nuclear charge of Z giving rise to the following non-relativistic Hamiltonian in atomic units {e = h = m = 47reo = 1) 

1 , 1 , Z Z 1 

^ = -:tV?--V + - (25) 

2 2 n J'2 ki-T-2| 

where r = \r\. With Z ~ 2 this is the Hamiltonian for the helium atom. Interestingly although the hydride atom {Z = 1) exists, 
and its energy may be calculated. Baker et al 1 1 1] showed that for nuclear charges below Z w 0.911 the system is unbound. 

Hooke's atom is an interacting system of two electrons held by a harmonic confining potential with Hamiltonian in atomic 
units of 

H = -^V^ - + + lu:\i + \ (26) 

2 2 2 2 \ri — r2\ 



Here lu is the characteristic frequency of the confining potential. The exact spatial entanglement and its approximations within 
LDA have previously been evaluated for Hooke's atom in Ref. 
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V. METHOD OF SOLUTION 



We solve the many -body Schrodinger equation, using 'exact' diagonalisation (EXD) with an efficient choice of basis. For two 
electrons in a spherically symmetric potential the ground state can only be a function of the distance of each electron from the 
origin and the angle between the electron vectors. Hence we may find the ground state using exact diagonalisation with a basis 
in three dimensions: 7-1, r2,0. 

We use the basis 



= R^{rl)RJ{r2)^J'^^Plicosi9)), (27) 
V47r V 47r 

where we employ the hydrogen-like wavefunction Ri{r) = Qi{r)e^°'^ . The choice of a is influenced by the strength of the 
nuclear charge Z. Here the Qi (r) are polynomials of degree i created via the Gram-Schmidt procedure such that the Ri are 
orthonormal, i.e., Ri{r)Rj{r)r'^dr = Sij. The Legendre polynomials Pi are orthogonal to each other and we explicitly 
include the normalisation factor for integration over all angles. We intend to calculate the matrix elements of the Hamiltonian 
with respect to the and (pi'j'i' basis functions Eq.|22] The external potential is spherically symmetric therefore the potential 
integrals become just ID integrals and are only non-zero when / = I'. This is propitious as they are the only integrals that have 
to be re-calculated each time in the iterative scheme. 
The kinetic term in the integrand is 

V?--— rf— (28) 



where Li is the angular momentum operator. As 

Pi{cos{e)) = Yi,n{eiAi)Y;^{02,<l>2), (29) 



m—-\--l 

(27Ti) 



then 

LlPi=l{l + l)Pi. (30) 

Hence the kinetic terms are also ID integrals in ri and r2 when /' = /, otherwise the integrals are zero. 

The Coulomb interaction integral is more complicated. First we expand the Coulomb interaction term using the Legendre 
polynomials 



fe=0 



where r< = min{ri, ^2} and ?■> = max{ri, r2}. 

Next we use a result Iil2il for an integration of three Legendre polynomials 

^ k 



Pi,iz)Pkiz)Piiz)dz^2l^Q Q . (32) 

If /' > k + I then the polynomial of degree k + I can be formed from Legendre polynomials each of degree less than I' 
hence the integral is zero, similarly for k > I + I' and I > k + I'. Here the 3j symbol (see for example Ref. [il3il ) is zero for 
I' + k + I = 2p + 1 where p is an integer. While for I' + k + I ~ 2p 



where 



(a + b-c)!(b + e-a)!(c + a-6)! 
Aiabc) = . (34) 

^ ^ (a + 6 + c+l)! 



Hence the integral for the Coulomb interaction when switching to spherical polar integration, including the normalisation, and 
integrating over angles is 
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i?j'(r2)i?j(r2)^dr2). (35) 



As the 3j symbol is zero for k > I + I' then there are at most / + I' double integrals, but the other constraints then produce 
a zero 3j symbol for many of these. Hence efficiency is increased: we have a small sum of double integrals rather than a triple 
integral to compute. 

We also consider the single particle KS equation (Eq.fTSll where fext(f") — —Z/r for helium-like atoms and Vf^{r) — 
for Hooke's atom. For approximating we use v^^^ = + We employ the fit of Perdew and Wang (al (see Eg. [23]) 
for v^;^^ while the expression for v]^^^ is that of Wigner and Seitz [f] (see Eq. \S2\ . We consider only spherically symmetric 
potentials so the KS equation reduces to a one dimensional problem. We solve this self consistently to obtain n'^'^*^*, the LDA 
approximation to the ground state electron density. 

This 'target' density is then used in the iterative scheme developed in Ref. fW] which allows to find the external potential for 
the interacting many-body system of ground state density n'^'^*", in this case the i-LDA system. This scheme is based on the 
relation 

<t.Hri) - ^-^|i?,:|N(ri) - n'-^^Xr-i)] + <,(ri). (36) 

Here Ei is the energy of the interacting system with external potential J^j ^ext(^j)' '^his many -body system is solved at 
every step by using the aforementioned method to give the many-body wavefunction from which we compute the den- 
sity rii. We aid convergence for both test-bed systems by mixing UgJ^ with 80% of w^^f, and iterate until the error 
/ d'^r |ni(r) — n'"S'^'(r)|/ J d^r n^^^^^^r) has reached a desired level on the chosen integration range. 

We note that the inversion scheme Eq.|36]could also be used to derive the unknown potential related to a generic density n'^^'^' 
obtained, e.g., from experimental measurements. 

VI. SPATIAL ENTANGLEMENT CALCULATION 

For the systems we are considering, the ground state wave-function can be factorised in orbital and spin part. The spin part 
is a singlet which is maximally entangled and constant. We therefore ignore the spin entanglement and investigate the spatial 
entanglement, i.e., the particle-particle entanglement stemming from the continuous spatial degrees of freedom. 

We use the two-electron wavefunction to calculate the linear entropy of the reduced density matrix 

L = 1 - Trp\ (37) 

as a measure of the spatial entanglement. Here pA — Tvb \tp) {tp\ is found by tracing out the system density matrix over the 
degrees of freedom of one of the two subsystems. The linear entropy was, perhaps, first put forward as a measure of purity of a 
quantum state by Zurek, Habib and Paz 1 14]. The purity of the reduced density matrix can be shown to give an indication of the 
number and spread of the terms in the Schmidt decomposition of the state and therefore a measure of the entanglement. Here, in 
the continuous case, we calculate the entanglement using 

PA(ri,r2) = J '^*{ri,ra)-^{r2,rs)dra, (38) 

Pi{ri,r2)^ J PA{ri,r3)pA{r3,r2)dr3, (39) 

and 

Trpl = J pl{r,r)dr. (40) 

We will investigate the change in spatial entanglement with the strength of the confining potential, which is characterised by 
uj or the nuclear charge Z for Hooke's and helium-like atoms respectively. 
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VII. HELIUM-LIKE ATOMS 

For helium-like atoms the (numerically) exact and the i-LDA entanglement are depicted in Fig. [T| There we see that the 
entanglement increases as the nuclear charge decreases which may be attributed to the increase in the ratio of the electron- 
electron interaction to the energy from the external potential. A similar situation was observed in the spatial entanglement 
of Hooke's atom fst]. We note that the spatial entanglement is low for the parameters considered, when compared with the 
theoretical maximum of L = 1. The 'exact' and i-LDA entanglement are close and the gap between them increases as the 
weight of the particle-particle interactions increase. This fits in with the LDA performing better when many-body interactions 
are weak: as the confining potential increases, many-body interactions will become less important relative to it, until the system 
may be fairly well approximated by a non-interacting Hamiltonian thereby giving rise to a product state solution. We note that 

0.025 



0.02 
0.015 

0.01 
0.005 



"1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

Z(e) 

FIG. 1 : The linear entropy L of the reduced density matrix versus the nuclear charge Z as a measure of entanglement for helium-like atoms 
and the corresponding i-LDA systems. 

although helium-like atoms have a bound state until Z w 0.911 we were unable to numerically solve the Kohn-Sham equations 
with the LDA satisfactorily for Z < 1.5. We hypothesise that for the exact solution with Z < 0.911 a transition to a resonance 
state — as studied recently for two electrons in a spherically symmetric square well ifisll and the spherical helium atom fT^ — 
would be expected to occur. This would result in a superposition of one electron constrained by the nucleus with the other 
unbound. Such a wavefunction would be expected to have a spatial entanglement of L = 1/2. 

VIII. HOOKE'S ATOM 

For Hooke's atom (Fig. |2]i we see that the accuracy of the entanglement when using the LDA is high, but decreases a little 
as ui decreases and interactions, as quantified by the ratio of the electron-electron interaction to the energy from the external 
potential, become larger. The value of the i-LDA entanglement for u = 0.00584, lower than the value for w — 0.00957, seems 
anomalous as does the change from an overestimate to an underestimate of the exact entanglement. We note that uj = 0.00584 
has the poorest match (0.053% error) between the i-LDA estimate and the actual LDA density of the uj displayed and is the 
smallest oj for which we could achieve a good match with a reasonable basis size of 8^^ wavefunctions. Therefore due to a 
seemingly high sensitivity of the entanglement to the wavefunction yet a relative insensitivity of the density to the wavefunction, 
we may have crossed a threshold so the small error in reproducing the LDA density is now not small enough to prevent a large 
error when calculating the LDA entanglement. However we also note that 0.053% error still represents a very good accuracy 
in reproducing the LDA density, and that when employing an even larger basis the difference between old and new density 
estimates was just 0.089%. A basis of 8"^ was also sufficient to calculate the exact entanglement to 3 significant figures. In 
addition we saw in previous work [3] that the difference between exact and LDA energies exhibited a non-trivial relationship 
with ljj, and a perturbed LDA-based approach to approximate Hooke's atom many-body wavefunction also produces a decrease 
of the entanglement in a similar range of ui (see Fig. O. These observations may imply that a non-trivial relationship between 
Lu and LDA-based estimates of the entanglement could be a real effect. We suggest that although the accuracy of the i-LDA 
entanglement for the two smallest values of uj may be questionable, the correct i-LDA entanglement for these values should not 
be too different and we are confident of the accuracy of the i-LDA entanglement for larger values of u. 
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FIG. 2: The linear entropy L of the reduced density matrix versus a; as a measure of entanglement for Hooke's atom and the corresponding 
i-LDA system. 

A. Comparison with other LDA-based entanglement estimates 

In previous work [3] we devised different ways to derive an LDA-based approximation to the entanglement. 

The first is a combination of gradient descent and evolutionary algorithms to derive an approximation for the i-LDA wave- 
function from which to calculate the entanglement. The main limitation of this approach is that the trial i-LDA wavefunction 
was based on a specific functional form. This variational-type approach included some parameters to be varied to improve 
performance. 

The second approach to derive an LDA-based approximation for the interacting many-body wavefunction is based on a 
perturbative scheme whose zero-order Hamiltonian includes the Kohn-Sham effective potential (and hence, through Vxc, partly 
the many-body interaction) with the LDA approximation for the exchange-correlation term. This approach improves as higher 
orders are taken into account, but a first-order wavefunction would still fail to reproduce the exact entanglement when many -body 
interactions become dominant. 

We see in Fig. [3] that, in the intermediate and strong particle-particle interaction regime, the exact i-LDA entanglement for 
Hooke's atom reproduces the system exact entanglement much more accurately than any of the approximations mentioned 
above. Approximations labelled as EAl and EA2 were achieved by using a combination of gradient descent and evolutionary 
algorithms [3J. In EAl we sought to minimise the difference between the density arising from the trial wavefunction and the 
exact density and then, within the wavefunctions giving an accurate enough match, we choose the wavefunction corresponding 
to the lowest ground state energy. The trial wavefunctions in this case are restricted to wavefunctions which can be separated 
in centre of mass and relative motion components. In EA2 we sought to minimise the difference between the density arising 
from the trial wavefunction and the exact density together with a minimisation of the energy of the trial wavefunction. The trial 
wavefunctions in this case used a finite number of functions in terms of ri, r2, and 6 from a general basis for two electrons 
in a spherically symmetric potential. The gradient descent allows the basis coefficients for the wavefunction that gives a local 
minimum to be found while the evolutionary algorithm allows more of the parameter space to be explored so in principle, given a 
large enough basis and enough time, the combination should find the global minimum. With our finite basis and limited run-time 
we see that, although this method gives the trend of the i-LDA entanglement, the approximate i-LDA entanglement from the 
evolutionary algorithm is much higher than the actual i-LDA entanglement. 

As mentioned above, in L3] we also considered a perturbative approach based on the LDA-KS equations, whose related 
interacting wavefunction converges toward the exact wavefunction as higher orders are considered. To first order we would hope 
for a good approximation to the exact entanglement that perhaps leans towards the i-LDA entanglement. We see in Fig. [3] that 
the entanglement results from first order becomes much greater than the i-LDA entanglement — but less so than the evolutionary 
algorithm — for stronger interactions, as quantified by the ratio of expectation of the interaction energy to the potential energy. 
However, the perturbation method allows an easier and more efficient way to approximate the entanglement within an LDA- 
related scheme, whose first order is accurate for low interactions; the exact i-LDA entanglement, although more accurate, is 
inaccessible without repeated complicated calculations: the inversion scheme requires solving the many-body equations more 
times than would be needed to ascertain the exact entanglement. 
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FIG. 3: The linear entropy L of the reduced density matrix as a measure of the entanglement obtained using the exact, i-LDA, first order 
perturbation of the LDA KS equations, and evolutionary algorithm (EAl and EA2) many-body wavefunctions. 



IX. POSITION SPACE INFORMATION ENTROPY 



We also investigate the position-space information entropy of the density 



Sn = — / n{r)\nn{r)dr. 



(41) 



as a proxy measure for the entanglement. This has been used to study the entanglement of the Moshinsky atom ITtIi and may also 
be thought of as a zeroth-order approximation to the von Neumann entropy of the reduced density matrix, where the off-diagonal 
terms of the reduced density matrix are set to zero [3]. As the von Neumann entropy of the reduced density matrix may be used 
as a measure of the entanglement and has been seen to behave similarly to L for Hooke's atom |3], it would be hoped that Sn 
should capture some of the behaviour of the entanglement. 



c 




FIG. 4: The position space information entropy Sn for the 'exact' and LDA densities for helium-like atoms. 

We see in Fig. |4] that Sn captures the monotonically decreasing trend of the exact and i-LDA entanglement. The Sn for the 
LDA is greater than the exact 5„ but the difference between the two does not seem to increase as Z decreases as is the case for 
the exact and i-LDA entanglement. The behaviour shown by 5„ is also more linear with increasing Z than the entanglement. 

For Hooke's atom we previously noted that Sn did not model the entanglement behaviour well over a much larger range of 
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confining potentials ' . We see in Fig.|5]that, over a smaller range of uj, Sn reproduces the general trend of the entanglement but 
Sn from LDA is always greater than the exact Sn which appears to not be the case for the entanglement (see Fig.|2]i. In addition 
we see that 5„ is almost a straight line when plotted against a; on a logarithmic scale, this is in contrast to the entanglement 
which can be seen to increase less fast as uj decreases even for the small range of w considered here. 
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FIG. 5: The position space information entropy Sn for the 'exact' and LDA densities for Hooke's atom. 



X. CONCLUSION 

We have calculated the spatial entanglement of the electrons in helium-like atoms corresponding to the exact system and to the 
use of the local-density approximation within density-functional theory using the concept of the interacting LDA system. The 
entanglement was seen to decrease with increasing nuclear charge Z. The maximum entanglement found was only a fraction ^ 
2.5% of the theoretical maximum. The entanglement of the i-LDA system was seen to reproduce the qualitative and quantitative 
behaviour of the exact entanglement with relatively good accuracy, although the discrepancy between the measures increased 
as the nuclear charge diminished. We also compared the exact spatial entanglement of Hooke's atom to the entanglement of 
the related i-LDA system and saw that the LDA was a reasonably good approximation to the entanglement, but the accuracy 
decreased as the strength of the confining potential decreased and the i-LDA entanglement became non-monotonic at very small 
confining potential. 

We also saw that previously considered approximations to the exact or i-LDA entanglement were not so reliable when the 
particle-particle interaction became dominant, but noted that they captured the trend of the entanglement, were successful at 
strong confining potentials and were often much simpler computationally. In particular the accuracy of the perturbative approach 
based on the LDA-KS equations could be systematically improved by going beyond first order. 

We considered the position space information entropy of the density, and show that it was able to efficiently give the general 
trend in the entanglement but not its detailed behaviour 

From our test-bed calculations on few-particle systems it would appear that the LDA may reproduce the entanglement be- 
haviour well, even for systems with a highly inhomogeneous density and fairly strong interactions. More studies are necessary 
to ascertain if this conclusion holds for many-particle systems. 



' Label on y-axis in f3], Figs. 5 and 12, should read 5^/ In 2. 
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